Regular and Anomalous Quantum Diffusion in the Fibonacci Kicked Rotator. 

G. Casati«, G. ManticaW and D.L. Shepelyansky^ 2 ) 
Wlnt. Ctr for the Study of Dynamical Systems, Universita dell'Insubria, Via Valleggio 11 1-22100 Como Italy 

and Istituto Nazionale Fisica della Materia, unita di Milano 

and Istituto Nazionale Fisica Nucleare, sez. di Milano 
< 2 ) Laboratoire de Physique Quantique, UMR 5626 du CNRS, 
Universite Paul Sabatier, F-31062 Toulouse Cedex 4, France 



o 
o 
o 



> 
m 

o 
o 
o 



c 

o 
o 



We study the dynamics of a quantum rotator kicked 
according to the almost-periodic Fibonacci sequence. A 
special numerical technique allows us to carry on this in- 
vestigation for as many as 10 12 kicks. It is shown that 
above a critical kick strength the excitation of the sys- 
tem is well described by regular diffusion, while below this 
border it becomes anomalous, and sub-diffusive. A law for 
the dependence of the exponent of anomalous sub-diffusion 
on the system parameters is established numerically. The 
analogy between these results and quantum diffusion in 
models of quasi-crystal and in the kicked Harper system 
is discussed. 



I. INTRODUCTION 

It can be reasonably argued that the quantum 
kicked rotator JjJ is the oldest, and more thoroughly 
investigated system in the realm of the so-called quan- 
tum chaos. Yet, we pretend in this paper that recourse 
to this system is still justified, and can provide new 
dynamical phenomena worth of investigation. Con- 
sider in fact the problem of transport in quantum ex- 
tended systems: for instance, a wave-packet is initially 
localized on a few states of a quantum lattice, and 
is then left free to evolve. A deep analysis carried 
out on model systems |^-|7]] with nearest neighbours 
interactions has shown that this dynamics is either 
localizedr-this is typically the case of random couplings, 
or, in the case of almost periodic couplings, anomalous 
(i.e. the diffusion coefficient is either zero or infinity- 
almost never finite) and intermittent [p|~[l4"||. 

It is well known that the kicked rotator for typical 
values of the parameters is akin to the first example 
just presented: its quantum dynamics is (eventually) 
localized in the action lattice. Therefore, it seems to 
feature a far simpler dynamics than the second class 
of systems mentioned above, or than its not-too-close 
relative, the quantum kicked Harper model [ pL5| , in 
which anomalous diffusion jl6j has been documented. 
A way to cure this simplicity is offered by an idea prof- 



itably employed in spin and few-levels systems 17 19| : 
rather than acting upon the system with a periodic 
impulsive force, we let the sequence of time-intervals 
between different "kicks" be almost-periodic. In few- 
level systems, this idea has lead to interesting spectral 
properties and dynamical behaviours. 

In applying this procedure to the kicked rotator we 
hope so to combine the complexity of an almost pe- 
riodic interval sequence and the richness of the infi- 
nite Hilbert space of the unperturbed problem, while 
maintaining at the same time a sufficient simplicity 
for the theoretical analysis: for this, we elect to in- 
vestigate the effect of a Fibonacci sequence of unitary 
evolution operators. 

This paper presents the results of a numerical anal- 
ysis of the system so defined. In previous studies of 
almost-periodic kicking of the quantum rotator, a lin- 
ear increase of the energy has been observed [ p0|j2l] , 
which amounts to regular diffusion. In alternative, 
power-law increase of energy with a non-unit expo- 
nent defines anomalous diffusion. We shall attach 
special importance to this exponent, which describes 
the growth of the number of the excited states of the 
system. Energy is just one of the indicators of the 
motion of this quantum system. We shall use here a 
larger set of indicators to describe, in the long time, 
the growth of the excitation of the system. Each of 
these quantities will be characterized by a power-law 
asymptotic behaviour, which in turn will define an as- 
sociated growth exponent. 

This analysis will show that two dynamical regions 
are present in this system. For large kicking strength 
the motion is diffusive: the growth exponents are 
equal to one half, and a proper diffusion coefficient can 
be defined. To the contrary, when the kick strength is 
less than a certain threshold, the motion typically fea- 
tures anomalous sub- diffusion: growth exponents still 
exist, but are smaller than one half. The investiga- 
tion is further carried on by varying the kick strength: 
in so doing, a power-law dependence of the growth 
exponents is observed. It is interesting to note that 
a similar situation has been found in quasi-crystals 
in two dimensions p2fl : in these systems, for strong 
hopping amplitude the wave-packet spreading over a 
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two-dimensional lattice is diffusive, while, in the op- 
posite case, spreading is characterized by anomalous 
sub-diffusion. 

This paper is organized as follows: in the next sec- 
tion we define the model. In Sect. Ill we introduce 
growth indicators and their power-law exponents to 
gauge anomalous diffusion in the dynamics of an ini- 
tially localized wave-packet. The numerical algorithm 
and the numerically computed long time dynamics of 
the system are presented in Sect. IV. Finally, we allow 
a random element in the definition of the dynamical 
system to verify the robustness of the obtained results, 
and to support the conjecture that anomalous diffu- 
sion in this model is induced by the almost-periodicity 
of the sequence of unitary operators. 



II. THE DYNAMICAL SYSTEM 

The basis of our investigation is the quantum kicked 
rotator this is the unitary evolution generated in 
the usual Hilbcrt space -^(S^tt) by the operators 



U(k,T) 



-ik cos 



(i) 



which depend on the two parameters k, T and where 
the momentum operator is n — —idg. The corre- 
sponding classical system is the well known Chirikov 
standard map |23|| , which is obtained letting k grow, 
and T go to zero, while keeping their product fixed. 
This product defines the classical chaos parameter, 
K = kT . In the classical case the system becomes 
globally diffusive for K > 0.9716.... 

Now, it is well known that for large values of K, and 
typical irrational values of classical and quantum 
dynamics of the kicked rotator are profoundly differ- 
ent: this is the content of the quantum phenomenon of 
dynamical localization, which has been well described 
in the literature, and which predicts that the quan- 
tum evolution of an initially focused wave-function is 
almost periodic in time, and localized on the momen- 
tum space of the unperturbed system at k = [|J . 

In this paper, we study a variation on this theme: 
while in the usual approach the unitary evolution of 
eq. ([I]) is acted repeatedly, at equally spaced times 
t n = nT, in our model the length of the time intervals 
in between kicks will not be constant. Equivalently, 
we can say that the operators U act now in Fibonacci 
sequence to produce the full quantum evolution U : 

U{t n ) = U(k, T n ) o U(k, T n -!) o . . . o U(k, n), (2) 

where t n — Y]"—-\ Tj is time, and where o denotes op- 
erator composition. As is immediately observed, op- 
erators in eq. (|^) are characterised by the same value 



of the quantum kick amplitude k, and differ solely by 
the value of T. The sequence of values of this latter 
is determined by letting Tj to take either one of two 
possible values, which for convenience we denote by A 
and B. We choose for A and B the values 

"X' B= *> 

where A = 1.3247... is an algebraic number, that solves 
the equation A 3 — A — 1 = 0. In a sense, A, B is the 
most irrational couple of numbers. 

We line up A's and -B's according to the well known 
Fibonacci sequence: this is obtained by an accre- 
tion rule on the initial finite pieces of the infinite 
sequence, o\ — A and oi — A, B. According to 
this rule, the next finite piece is obtained by joining 
these two sequences, the first at the end of the second: 
CT3 = A, B, A. In general, 



OVi+l — &n U CT n _l, 



(3) 



where the symbol U means here joining the right and 
left sequences in the specified order, so that 04 = 
A, B, A, A, B, and so on. It is then easy to see that 
the length of o~ n is equal to the n-th Fibonacci num- 
ber, F n . Eq. (||) expresses a property that will be 
instrumental in the following to obtain the quantum 
evolution up to large times. 



III. GROWTH EXPONENTS OF AN 
INITIALLY LOCALIZED WAVE-FUNCTION 

We now let ip(0) be the initial state of the evo- 
lution, and we compute the time-dynamics ii(t n ) = 
U(t n )ip(0). Typically, we choose for initial state ip(0) 
the usual eo momentum basis state of -^(SW), the 
full basis being e n = e m6 , with n E Z. 

As we have already remarked, the quantum kicked 
rotator, for typical irrational values of the kicking pe- 
riod T, displays the phenomenon of quantum localiza- 
tion: action in the system eventually saturates in time, 
and quantum diffusion stops. On the other hand, for 
kick sequences with random time intervals after each 
kick, the increase in momentum is diffusive and never 
Of particular interest is then 
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to study the kicking sequence defined in the previ- 
ous section, which is neither periodic, nor random: 
it belongs to the family of almost-periodic sequences. 
We shall see that under these circumstances dynami- 
cal localization leaves places to a derealization which 
shares many characteristics with that taking place in 
almost-periodic lattice systems. 

To study the dynamics of the quantum wave- 
function, it is convenient to introduce a number of 
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quantities suitable to define its spreading: we shall 
call them collectively growth indicators. The first ob- 
vious choice are moments. Let v a (t) be the moment 
of index a of the probability distribution over the lat- 
tice {e„} given by |(e„, tp(t))\ 2 , where obviously (•,•) 
denotes the scalar product in ^(SW): 

v a {t):=Y,\{e n ,m)\>\ a - (4) 

n 

For instance, v 2 is the usual second moment-the 
dispersion- which, in the case of regular diffusion, 
grows linearly in time. In general, we cannot expect 
u 2 to behave as t, but super-diffusive and sub-diffusive 
growth will be the rule. For this, we also define the 
growth exponent function (3 by the asymptotic relation 

v a {t) ~ t a(s(a) for t^oo, (5) 

where the asymptotic behaviour is to be understood 
in a suitable sense. The function [3(a) is the so-called 
quantum intermittency function, described at length 

in m. 

When a tends to zero, the moment v a tends to the 
constant uq = 1. Nonetheless, one can define the limit 
of the function (3 for a tending to zero, /3(0), a number 
which also appears in the scaling of the logarithmic 
moment: 

^|( en! ^))| 2 log(|n|)~/3(0)log(t). (6) 

In the quantities (|^) and (^| the ordering of the ba- 
sis is relevant: a function of n appears in the summa- 
tions. We can nonetheless define and compute quan- 
tities which do not make reference to this characteris- 
tics. The first of these is the entropy: 

-^|(e n ,^(i))| 2 log(|(e, i ,V(t))| 2 )^7 + /3elog(i), 

n 

(7) 

whose asymptotic behaviour (here and in the following 
to be understood in the infinite time limit) defines 
the exponent [3 e . In the case of localization, (3 e is 
null, and e 7 (the constant contribution in eq. (]?])) 
gives a measure of the localization length. Another 
basis independent quantity is the inverse participation 
ratio, 

J2\(e n ,m)\ 4 -t-^, (8) 

n 

with its exponent [3 p . 

We notice that, in principle, the exponents [3 of all 
the quantities so far defined are different. Yet, when 



finite, they are all consistent with the fact that the 
spreading of the wave-packet over the basis takes place 
in a power-law fashion. This is indeed the case: in fig. 
(|l|) we plot the ^ power of a set of moments v a versus 
time, in doubly logarithmic scale, together with one 
over the inverse participation ratio, and with the expo- 
nential of the entropy and of the logarithmic moment. 
It is immediately appreciated that quantum diffusion 
in this case is regular: these curves are parallel to the 
reference line \ft, so that all growth exponents are 
equal to one half. 

To render even more evident the fact that we are 
in the presence of regular diffusion we plot in fig. (|2|) 
a snapshot of the wave-packet at a large time, for the 
same case of fig. (0): The fitting line is a normal 
distribution. 



IV. DETECTING ANOMALOUS DIFFUSION 
IN LARGE TIME DYNAMICS 

Let us now investigate the effect of varying the kick 
strength k: in particular, let us decrease it from the 
value k = 10 at which we observed regular diffu- 
sion. Two difficulties arise when performing this anal- 
ysis: since we observe a reduced derealization of the 
wave-packet, very large time-scale computations are 
required to say something reliable about the asymp- 
totic behaviour. Secondly, large systematic oscilla- 
tions in the behaviour of the moments (and of the 
other quantities) tend to overwhelm and submerge the 
power-law increase, the more the slower this increase. 

Let us first present the cure to the second difficulty: 
we can introduce a magnetic flux a in the evolution 
(|]J), which is nothing else than replacing the operator 
h with n + cf>: is a real number, so that for = 
we obtain the usual kicked system. By performing 
the evolution with a set of values of <fi, and by aver- 
aging the results with respect to this sample, we can 
reduce the systematic fluctuations super-imposed to 
the leading power-law behaviour. Please observe that 
the classical dynamics is left invariant by the intro- 
duction of this flux. 

For the first difficulty, a more complex strategy is 
required. We can push the evolution to extremely 
large times by resorting to a matrix algorithm. Let 
U Fn be the unitary evolution operator acting from 
time zero up to the time of the kick number F n , the 
n-th Fibonacci number. Then, eq. (^|) implies that 
these operators are simply related by 

where o is the usual operator multiply, and where the 
order of operators is essential, as it was in eq. (0). 
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This relation must be initialised by setting 

U Fl = U(k,A), 

U F2 = U(k,B)oU(k,A), 



(10) 



so to produce the full sequence {U F ™}. In fact, the 
matrix elements of U (K, r) are explicitly known, and 
eq. (j^) involves then only a matrix-to-matrix multipli- 
cation. Of course, the cost of this operation scales as 
the cube of the size of the matrix, and this is certainly 
more expensive than the usual evolution effected via a 
sequence of fast Fourier transforms. Nonetheless, this 
technique becomes advantageous when dealing with 
extremely large times, for in n steps one reaches up to 
F n map iterations, and it is well known that Fibonacci 
numbers are geometrically increasing in n. 

We therefore effect the matrix multiplication for 
n = 1, . . . , N and then look at the zeroth coloumn 
of the resulting operator U Fn : this is the evoluted of 
the initial state eo- In our numerical experiments, the 
initial state of the evolution will always be eo- Suitable 
combinations of the columns can provide the evoluted 
of any arbitrary initial state, if desired. 

A question of numerical concern is to be discussed 
prior to showing the results of these calculations: 
when dealing with equations (|8|Jl0|) a finite truncation 
of the operators involved is necessary. Now, for any 
practical purpose U(k, r) is a banded matrix, and we 
must obviously choose the size of the numerical trun- 
cation much larger than the band size. Yet this is not 
enough: appropriate boundary conditions at the ba- 
sis edges must be imposed. Two choices are at hand: 
Dirichlet, and Neumann. In the latter case, we are 
producing the unitary evolution on the torus, while in 
the first case we end up with a non-unitary evolution. 

Having in mind what we want to obtain (i.e. the 
evolution of the eo basis state in the full cylinder 
space) we have adopted Dirichlet conditions, for two 
reasons. First, when the truncation size is suffi- 
ciently large, both Dirichlet and Neumann conditions 
must produce an excellent approximation of the true 
U Fn eo- Secondly, we can gauge numerically two quan- 
tities: the normalization of the zeroth column of 
U Fn (which is a common technique) and its effective 
dimension-which is less common, but more instruc- 
tive. In fact, this latter was defined above as the num- 
ber of components of U Fn eo of magnitude larger than 
a certain threshold: this quantifies the dimension of 
the wave-packet after kicks, and can be profitably 
used to control the optimal size of the truncation. It 
must be added that the effective dimension in itself 
has physical and mathematical relevance 

In summary, we have been able to reach easily about 
sixty iterations of the Fibonacci multiplication, eq. 
(|l0|), which correspond to more than 10 12 usual iter- 



ations, and this with a basis size of the order of the 
thousand. Numerical stability, controlled by various 
techniques, is here the limiting factor. 

A set of typical results is shown in fig. (j^), where we 
plot the second moment, v^, as a function of time, in 
doubly logarithmic scale, for a set of different values of 
k. We clearly see that these values decrease with the 
latter, but, more important, the slopes of the curves 
(i.e. the growth exponent (3(2)) also decrease with k. 
As a consequence, we can safely conclude that quan- 
tum diffusion is anomalous in this model for a wide 
set of parameters S . 

At this point, a word must be spent about intcrmit- 
tency: this is present when the function (3(a) (see eq. 
(H)) is not constant. In the case at hand, the variation 
of (3 with a is smaller than the uncertainty with which 
(3 is determined, both due to numerical effects, and 
to the super- imposed oscillations of the moments v a . 
We can therefore only conclude that intermittency, if 
present, is low. 

This fact is not totally negative: by averaging over 
a, we can define a more reliable growth exponent (3 av 
which is now function only of the kick amplitude k. 
This immediately prompts for the study of this de- 
pendence. Fig. (ji|) plots the numerical results ob- 
tained by the procedure exposed above. Two regions 
clearly emerge from the investigation of this picture: 
for large values of k we observe the regular diffusive 
value (3 av — \- Anomalous diffusion is observed for 
smaller values of k: quite interestingly, we find a 
power-law behaviour of the growth exponent of the 
form 



(3av(k) ~ k , 



(11) 



with t] very close to the value |. The transition from 
anomalous to regular diffusion takes place for k > k c : 
for the set of values we have chosen, this critical value 
k c is approximately equal to two. 

Let us now investigate the origin of anomalous dif- 
fusion in this model. When acting on the basis set 
e n — e m6 the free evolution e l ~ 9 » produces the phase 
factor e" 1 !™ . The arithmetic nature of T/2n is at 
the root of the spectral properties of the usual kicked 
rotator. The numerical studies presented so far have 
been carried on for the most irrational pair A, B. We 
have found similar results for other irrationally related 
pairs. We are so lead to conclude that the nature of 
the observed unbounded diffusion lies in the almost- 
periodic arrangement of in-kicks intervals A and B. 

To substantiate this hypothesis we have replaced 
the phases -jn 2 by two sequences of equally dis- 
tributed pseudo-random numbers, one for the oper- 
ator U(k,A) and one for its companion U(k 1 B). In 
so doing, we have obtained quite similar results as 
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those reported above. In fig. || the exponents f3 av for 
this case are also reported. 



V. CONCLUSIONS 

We have studied the dynamics of a quantum rota- 
tor kicked at times generated by the almost periodic 
Fibonacci sequence. Contrary to that of the usual 
kicked rotator, this dynamics nevers shows quantum 
localization. We have introduced and computed var- 
ious indicators of the spreading of an initially local- 
ized wave packet, which have permitted to show that 
the dynamics features regular diffusion for large val- 
ues of the kick amplitude (k > k c ) and anomalous 
sub-diffusion for small values (k < k c ). In this latter 
range, the exponent of this diffusion appears to fol- 
low a power-law behaviour with the kick amplitude, 
Pav(k) ~ k~s . This relation, although only numeri- 
cally established, is quite interesting, and deserves in 
our opinion further investigation. 

Similarities and differences between the Fibonacci 
kicked rotator and other quasi-periodic models are to 
be noted. We have already observed that anomalous 
diffusion is typically found in almost periodic, one- 
dimensional lattice systems. It also appears in the 
kicked Harper model. To the contrary, the Fibonacci 
kicked rotator is qualitatively different from the rota- 
tor acted upon by equally spaced kicks, with an ampli- 
tude k which is a quasi-periodic function of time. The 
case when this function contains m incommensurate 
frequencies has been studied in [f26|-p8| . It was shown 
that this model is a dynamical analog of the Anderson 
localization in a space of effective dimension d = m+1. 
In this way the usual kicked rotator (m = 0) corre- 
sponds to d = 1, and-of course-is always localized. 
For to = 1 the excitation is still always localized, but 
the localization length grows exponentially with k 26 1 , 
in analogy with the Anderson localization in d = 2. 
Finally, for to > 1 i.e. d > 2, a transition from local- 
ization to diffusive excitation takes place above some 
critical kick amplitude, in analogy with the Anderson 
transition for d > 2 §7|J||. As we see this behaviour 
drastically differs from the one we find in this paper 
for the Fibonacci kicked rotator. 

In our opinion this difference might be due to the 
fact that in [p6|-p8| the kick amplitude k is an analytic 
function of incommensurate phases (frequencies). For 
this, such a model can be mapped into an effective 
solid-state hopping model, similar to the Anderson 
model, with hopping only to a finite number of nearby 
sites. On the contrary, the almost-periodic sequence 
of unitary operators of the Fibonacci kicked rotator 
renders the situation much richer, and gives rise to a 



transition from regular to anomalous diffusion. 

As a matter of facts, the dynamics of the Fibonacci 
kicked rotator seems more similar to a wave spread- 
ing on a two-dimensional quasi-crystal lattice. In- 
deed, studies of quantum diffusion over an octagonal 
quasipcriodic tiling have shown a similar transition 
from anomalous to regular diffusion p2|] . However, in 
spite of this first glance similarity more detailed stud- 
ies are required to establish a quantitative relation be- 
tween these models, and to gain a better theoretical 
understanding of the results presented in this paper. 
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FIG. 1. Regular quantum diffusion. Growth indicators 
are plotted versus time given in number of kicks (in ab- 
scissa) for the case k = 10. The lowest curve (dots-dashes) 
is the reference line yfi. The other curves are (from bottom 
to top): the exponential of the entropy (long dashes), the 
square root of the second moment (short dashes) , the sixth 
root of the sixth moment (dots) , one over the inverse par- 
ticipation ratio (dots-dashes), and the exponential of the 
logarithmic moment (continuous line). 




FIG. 2. Wave packet for the case of fig. 1. Plotted 
are occupation probabilities p n — \(e n ,ip(t))\ versus site 
number n at time t — 1.797 10 4 (thin line). The thick line 
is an interpolating normal distribution. 
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FIG. 3. Second moment u 2 as a function of time, 
for k = 0.125,0.25,0.5,1,2,4. (bottom to top diamond 
curves). Also shown are the fitting lines from which the 
exponent /3(2) is extracted. 
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FIG. 4. Growth exponent /3 a v (crosses) obtained as the 
average of f3(a), for a = 2, 4, 6 , versus the kick amplitude 
k. The horizontal line marks the value f3 av — | ; the fitting 
line for small values of k grows proportional to k v , with 
7] = g. Also reported (diamonds) are the result obtained 
when random rotation phases replace the irrational phases 
Tn 2 /2 (see text). 



